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ABSTRACT 

We present a simple approach for obtaining robust values of astrophysical parameters 
from the observed colour-magnitude diagrams (CMDs) of star clusters. The basic 
inputs are the Hess diagram built with the photometric measurements of a star cluster 
and a set of isochrones covering wide ranges of age and metallicity. In short, each 
isochrone is shifted in apparent distance modulus and colour excess until it crosses over 
the maximum possible Hess density. Repeating this step for all available isochrones 
leads to the construction of the solution map, in which the optimum values of age and 
metallicity - as well as foreground/background reddening and distance from the Sun 
- can be searched for. Controlled tests with simulated CMDs show that the approach 
is efficient in recovering the input values. We apply the approach to the open clusters 
M 67, NGC 6791, and NGC 2635, which are characterised by different ages, metallicities 
and distances from the Sun. 

Key words: (stars:) Hertzsprung-Russell and colour-magnitude diagrams; (Galaxy:) 
open clusters and associations: general. 


1 INTRODUCTION 

Colour-magnitude diagrams (CMDs) encapsulate important 
properties of single or composite stellar populations. The 
most relevant are the age, metallicity, binary fraction and, 
perhaps, the initial mass function. These are followed by 
the foreground/background reddening, distance from the 
Sun, differential reddening (especially for embedded clus¬ 
ters), and even age spread (for young clusters) of the con¬ 
stituent stars. Combined, these parameters contribute - in 
different degrees - to shaping and defining the CMD mor¬ 
phology of a star cluster. Obviously, the definition of the 
evolutionary sequences of a single stellar population in a 
CMD increases with the number of member stars. 

A proper or robust determination of astrophysical pa¬ 
rameters of star clusters is an important source of con¬ 
straints for tracing the Galactic structure and understanding 
its dynamics. Fundamental parameters of young star clusters 
are important also for studies of dynamical state and cluster 
dissolution time scales (e.g. Goodwin||2009 Lamers^_Baurtt 
gardt fc Gieles||2010|, infa nt mortalit^(J{e.g. Lada & Lada 
2003 Goodwin &: Bastian||2006 1, star-formation rate in the 


Galaxy (e.g. Lamers fc Gieles||2006 Bonatto fc Bica||2011 |, 

among others. This kind of analysis started with Trumpler 
||1930[ ) and continues to date (e.g. |Gamargo, Bica fc Bonatto 
2013 and references therein). 


In this context, the existence of methods to help in 
the extraction of parameters from CMDs is fundamental. 
The simplest approach is to compare a given CMD with 
those built with stellar populations of previously known pa¬ 
rameters. Depending on the quality of the CMD (basically, 
the definition of the evolutionary sequences), this exercise 
should yield relative values of the age, distance from the 
Sun and metallicity. A more quantitative approach is trying 
to find one isochrone (i.e., a population of stars having the 
same age and metallicity but different masses) that fits the 
overall CMD morphology of a star cluster. In principle, this 
can be - and has been - done by eye, but usually restricted to 
a small number of isochrones. This method is not practical 
for large sets of isochrones. 

However, a deeper interpretation of CMDs in terms of 
intrinsic cluster properties such as mass, age, metallicity, 
binary fractions, etc, requires more than estimates or rela¬ 
tive quantities. Thus, several approaches have been built for 
extracting reliable parameters from CMDs of star clusters. 
For instance, a Bayesian technique to invert CMDs of main- 
sequence and white dwarf stars has been intr oduced by|von 
Hippel et al.| ( |2006[) and furth er dev eloped by|DeGennaro et 
al. 120091 and van Dyk et al. (20091. Naylor &: Jeffries j2006 1 
apply maximum-likelihood statistics to Hess diagranjj sim- 


^ Essentially, the early (fO ~ 40Myr) phase in cluster evolution 
affected by significant disruption processes. 


^ Hess diagrams ( |Hess| |1924| contain the relative density of 
occurrence of stars in different colour-magnitude cells of the 
Hertzsprung-Russell diagram. 
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ulations (including binaries) to derive distances from the 
Sun and ages. However, they do not consider differential 
reddening, and their method appears to be more efficient 
for clusters older than ~ 30 Myr. |Hillenbrand, Bauermeister] 


& White (20081 model CMDs of star-forming regions and 


young star clusters by means of varying star-formation his¬ 
tories. da Rio, Gouliermis & Gennaro (20101 include differ¬ 
ential reddening, age spreads, and pre-main sequence stars 
to the Naylor & Jeffries (20061 method, but have to adopt 
distance from the Sun and reddening values from previous 
works. Kerber et al. (20021 use likelihood statistics to iden¬ 
tify the synthetic CMD that best reproduces the observed 
one, with membership probabilities determined from com¬ 
parison among CMDs extracted from the cluster and a con¬ 
trol field. More recently, Bonatto, Lima & Bica (20121 sim¬ 
ulated CMDs including cluster (stellar) mass, age, distance 
modulus, star-formation spread, foreground and differential 
reddening, and binary fraction to reproduce the observed 
CMDs of young star clusters. 

In general, the above approaches are designed to provide 
a thorough description of the stellar population properties 
giving rise to an observed CMD. Unfortunately, they usually 
invoke many free parameters and assumptions; besides, they 
also tend to overlook some degree of degeneracy associated 
with, for instance, differential reddening, age spread and bi¬ 
naries. They also do not account for the fact that the low 
photometric completeness affecting mainly the faint (low- 
mass) stars - especially at the central region of populous 
star clusters - may mimic a flat mass function, thus leading 
to confusion between actual dynamical evolution and com¬ 
pleteness effects. 

In many - if not most - cases, the interest is just to ex¬ 
tract reliable values for the age, metallicity, reddening and 
distance from the Sun of a star cluster, based only on the ob¬ 
served CMD and a set of isochrones. That’s the perspective 
we’ll explore in this work. 

This paper is organised as follows: in Sect. we de¬ 
scribe the new isochrone fitting approach. In Sect.I^we test 
the parameter recovery efficiency with simulated CMDs. In 
Sect. 1^ we apply the approach to actual star clusters and 
discuss the results. Concluding remarks are given in Sect.[^ 


2 FITTING THE HESS DIAGRAM 

The first step of the approach consists in building the Hess 
diagram from the photometric measurements of the ob¬ 
served stars in a given star cluster. As a compromise between 
resolution and computational time, the Hess diagrams used 
in this work are built with magnitude and colour cells of di¬ 
mension A mag — 0.02 and A col = 0.01, respectively. The 
smearing effect introduced by photometric uncertainties is 
explicitly considered in the Hess diagrams. Assuming that 
uncertainties follow a normal distribution around the mean 
values, we compute the magnitude and colour fraction (sum¬ 
ming for all stars) that occurs in a given cell of the Hess di¬ 
agram. In practice, this value corresponds to the difference 
of the error functions computed at the cell borders. Per¬ 
haps, this rather continuous nature is the most important 
advantage of Hess diagrams over the discreteness of CMDs. 
By definition, the sum of the colour and magnitude density 
over all Hess cells is the number of input stars. 



Figure 1. Each point represents an available isochrone built with 
the corresponding metallicity and age. 


In short, the present approach consists in finding the 
isochrone (or range of isochrones) that best maps the high- 
density ridge occurring along the evolutionary sequences. 
Quantitatively, this refers to the isochrone that crosses 
through the largest possible number of high-density cells in 
the Hess diagram. This measure is quantified by means of the 
fit index (find), which contains the total Hess-density sum 
of a given isochrone. Specifically, for an isochrone (charac¬ 
terized by a given age and metallicity) displaced on the Hess 
diagram by some amount of distance modulus and colour- 
excess, the fit index is quantified as 


find — ^ ^ H{m,i , Ci ), 


( 1 ) 


where N is the number of magnitude points in the isochrone, 
mi and a are the magnitude and colour corresponding to 
the ith isochrone point, and Hirrii, a) is the respective Hess 
density. Additionally, this procedure naturally leads to lower 
values of find for secondary sequences (i.e. lower-Hess den¬ 
sity) occurring in, e.g. high binary fractions. 

Given the wide-range coverage in age, metallicity, and 
mass - and the relative simplicity in obtaining them - we 


work here with the PARSEC isochrone set (Bressan et al. 
20121, version 1.2^ These isochrones are computed for 
a scaled-solar composition and follow the relation Y = 
0.2485 -I- 1.78Z. They consider the solar metal content as 


® available at http://stev.oapd.inaf.it/cgi-bin/cmd 
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Zq = 0.0152. It should be noted that our approach works 
with any isochrone set, as long as the isochrones cover the 
ranges of age, metallicity and stellar mass expected to occur 
in star clusters. 


Regarding ages and the total metal content Z, we con¬ 
sider isochrones distributed within 0.0001 ^ Z ^ 0.03 and 
ages in the range 1 Myr ^ tA ^ 13 Gyr, with higher reso¬ 
lution both for lower metallicity and younger ages. In sum¬ 
mary, we work here with a base of 2808 isochrones. For clar¬ 
ity, the distribution of isochrones on the Zx Age plane is also 
shown in Fig.[^ Obviously, the approach can be equally used 
with the full or restricted ranges of metallicity and age. 


Once the metallicity and age ranges are selected, each 
isochrone belonging to those ranges is used to search for the 
values of apparent distance modulus (DM) and colour ex¬ 
cess (CE) that produce the best fit to the Hess diagram. In 
the present context, a ht represents the density sum over all 
Hess cells crossed by the isochrone (see eq. [^. The search 
ranges for DM and CE must also be given at the start. Thus, 
our procedure reduces to two the dimensionality of the prob¬ 
lem, i.e., finding the best values of DM and CE. This might 
lead to losing some information on the error propagation in 
both parameters. However, as discussed in Sect. the final 
errors in DM and CE (and the remaining parameters) are 
computed based on the morphology of the plane ZxAge as 
a function of find- The search for the best values of DM and 
CE is carried out by the global optimisation method known 
as adaptive simulated annealing (ASA), which is relatively 
time efficient and robust (e.g. Goffe, Ferrier fc Rogers|1994 1. 
ASA is able to distinguish and escape from different local 
maxima in its search for the absolute maximum (e.g. |Bon-| 
atto, Lima fc Bica|2012|. 


As a technical caveat, we remark that an isochrone is 
a discrete set of stellar mass, magnitudes and colour values. 
But, by construction, the number - and consequently, distri¬ 
bution - of these points is not the same for all isochrones. So, 
isochrones having a higher number of points would naturally 
end up being favoured by the present approach. To avoid 
this, each isochrone is regularised at the start, so that the 
line connecting 2 consecutive isochrone points is discretised 
with the same spacing as that adopted for the Hess cells; 
additionally, each cell crossed by an isochrone is counted a 
single time. 


The end result of applying the above procedure to all 
isochrones (belonging to the selected age and metallicity 
range) is the solution map, which contains the fit index for 
all values of Z and age. For visualisation purposes, this map 
is first normalised by the find corresponding to the best so¬ 
lution. Then, the map resolution is increased by means of a 
bi-cubic spline interpolation. The latter step is necessary to 
find the signihcant solutions and compute their uncertain¬ 
ties (Sect.[^. In the context of the present work, a solution 
would correspond to a region (hot-spot) in the map where 
the highest values of find are concentrated. Then, the opti¬ 
mum values of age, metallicity, apparent distance modulus, 
and colour excess (and, consequently, distance from the Sun) 
can be derived from the hot-spot morphology. 



Figure 2. Top: ZxAge solution maps for Models#! and #2. The 
fit index scale is shown by the colour box (top-right). Bottom: V X 
(B — V) Hess diagrams together with the best-solution isochrones. 
Isochrone legends show the metallicity (Z), age (Gyr) and distance 
from the Sun (kpc). 


3 TESTING ON SIMULATIONS 


In this section we employ simulated CMDs to investigate the 
efficiency of our approach in retrieving the input values of 
relevant astrophysical parameters under very different num¬ 
bers of stars. The model CMDs correspond to relatively old 
star clusters with total stellar masses of 250 Mq (Model#l) 
and 10000 M0(Model#2), located at 1.1 kpc from the Sun. 
For the sake of realism, we assume that at this distance, 
B and V photometry obtained with typical ground-based 
telescopes would reach apparent magnitudes ~ 20 still with 
relatively small errors. 

For consistency with the remainder of this work, pho¬ 
tometric errors at given B and V magnitudes are simulated 
following the average error distribution as a function of mag¬ 
nitude observed for the star clusters M67, NGC6791 and 
NGC2635 (Sect. [^. Formally, for stars fainter than B = 10 
and V = 10, the corresponding photometric uncertainties 
are computed according to as = 0.01 -I- 1.0 x 
and av = 0.01 + 1.5 x ; brighter stars have 

(Js = 0.01 and ay ~ 0.01. The simulated stellar mass dis¬ 


tribution follows a Kroupa (20011 mass function, with each 


star being attributed B and V magnitude values according 
to its mass. The mass/light relation was extracted from a 
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specific PARSEC isochrone built according to the Bessell 
& Brett (Bessell 1990 Bessell & Brett 19981 photometric 
system. 


Except for the total stellar mass, the models are identi¬ 
cally built with the isochrone corresponding to 3.5 Gyrof age 
and total metallicity Z — 0.012; additional model parame¬ 
ters are the colour excess E{B — V)= 0.21, and the binary 
fraction /i,i„ = 0.3. At 3.5 Gyr, the corresponding isochrone 
contains stars in the mass range 0.09 ^ m{MQ) ^ 1.4. For 
more realism, photometric errors and scatter are also added 
to the stars. Models^l and ^2 end up having ~ 100 and 
~ 5000 stars brighter than V = 20, respectively. The differ¬ 
ence in number of stars, as well as the restricted photometric 
range, are useful for testing the efficiency of the approach 
under different morphology definition of the evolutionary se¬ 
quences on the corresponding Hess diagrams. 


To test the approach we considered all isochrones hav¬ 
ing age and metallicity within the ranges 100 Myr — 13 gy 
and 0.001 — 0.030, respectively. Fig.(top panels) shows the 
normalised, high-resolution solution maps for the three mod¬ 
els. Glearly, the convergence pattern on the age x metallicity 
plane gets tighter as the number of stars (and so, the def¬ 
inition of the evolutionary sequences on the Hess diagram) 
increases, which should be reflected both on the number of 
secondary solutions and the scale of the fit parameter un¬ 
certainties (see below). As expected, the best solution corre¬ 
sponds to the point with the highest find in the map. After 
locating it, we fit the map with a two-dimensional Gaus¬ 
sian, with axes along the Z and age directions; the values 
of age and Z at the maximum find are kept fixed. The free 
parameters, the Gaussian amplitude together with Z and 
age dispersions (ctz and cr/i) are searched again with ASA. 
The best-fit values of az and cta are taken as the global 
uncertainties of Z and age. This Gaussian is subsequently 
subtracted from the map, and the new (if any) maximum 
is searched for, thus giving rise to a set of solutions ranked 
by the respective find- Each solution corresponds to a single 
value for DM and CE. Then, considering all solutions occur¬ 
ring within ±(Tz and ±(Ja, we compute the weighted (using 
find as weight) average and Icr dispersion for both DM and 
CE. As a caveat we remark that, although this approach has 
the advantage of allowing the detection and characterization 
of multiple solutions, the find morphology around the solu¬ 
tions in Figs.[^and[^clearly deviate from a simple Gaussian- 
ity. The extent to which the deviation affects the parameter 
uncertainties and alternative approaches for finding errors 
will be explored in forthcoming work. The simulated CMDs 
and respective Hess diagrams are also shown in Fig. [^(mid¬ 
dle and bottom panels), together with the isochrones and 
parameters corresponding to the best solutions. Since the 
optimum Z and age are searched for on the interpolated, 
high-resolution map, the numerical values returned do not 
necessarily coincide with those of the input isochrones. 


The best-fit parameters obtained in the simulations 
are given in Table where we restrict to solutions hav¬ 
ing find > 0.7. Reflecting the low number of stars and, 
thus, less constrained CMD (and Hess diagram) morphol¬ 
ogy, Model^l has 3 solutions with find > 0.7. But, they all 
are consistent with the input values, within the uncertain¬ 
ties. 


In summary, tests with simulated CMDs show that our 
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Figure 3. Analysis of M67 (top panels) and respective model 
(bottom). Age and metallicity values may differ from those in 
Table [^because they correspond to those of the nearest-available 
isochrones. The fit index scale is shown by the colour boxes 
(right). Isochrone legends show the metallicity (Z), age (Gyr) and 
distance from the Sun (kpc). 


approach is efficient in recovering the input parameters and 
providing reasonable uncertainties. 


4 TESTING ON STAR CLUSTERS 

At this point we test our approach with actual star clus¬ 
ters. For this exercise we selected the relatively old open 
clusters M67 (NGC2682) and NGC6791, and the younger 
NGC 2635. They were selected for having relatively recent 
and uniform CCD photometry easily accessible through 
VizieFij^ In addition, they are characterised by different 
ages, metallicities, number of stars, distances from the Sun 
and amount of field-star contamination. Both M 67 and 
NGC 6791 are older and closer to the Sun and, perhaps, 
more populous than NGC 2635. One consequence is that 
both older clusters produce clearly discernible evolutionary 
sequences in CMDS or Hess diagrams (Fig.[^. NGC 2635, on 
the other hand, has a more ambiguous (less populated and 
more contaminated) CMD morphology. Thus, these three 


http://vizier.u-strasbg.fr/viz-bin/VizieR 
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Table 1. Parameters for models and star clusters 


find 

(1) 

Age 

(Myr) 

(2) 

Z 

(3) 

E{B-V) 

(mag) 

(4) 

d© 

(Kpc) 

(5) 

Model#l (109 stars) 

1.00 

3462 ± 160 

0.010 + 0.001 

0.23 + 0.01 

1.07 + 0.02 

0.85 

3173 ± 290 

0.009 + 0.001 

0.26 + 0.02 

1.07 + 0.03 

0.83 

3519+ 135 

0.013 + 0.001 

0.20 + 0.02 

1.11 + 0.02 

Model7^2 (5011 stars) 

1.00 

3500+ 120 

0.012 + 0.001 

0.22 + 0.01 

1.08 + 0.02 

M 67 (446 stars) 

1.00 

3292 + 300 

0.008 + 0.002 

0.16 + 0.03 

0.79 + 0.03 

0.82 

3542 + 150 

0.019 + 0.003 

0.02 + 0.02 

0.88 + 0.02 

NGC 6791 (2092 stars) 

1.00 

7694 + 505 

0.019 + 0.001 

0.23 + 0.02 

3.87 + 0.07 

0.97 

7472 + 633 

0.028 + 0.002 

0.17 + 0.01 

4.03 + 0.06 

0.85 

7556 + 250 

0.025 + 0.001 

0.18 + 0.02 

4.02 + 0.06 

0.78 

8000 + 440 

0.012 + 0.001 

0.29 + 0.02 

3.67 + 0.07 

NGC 2635 (481 stars) 

1.00 

500 + 100 

0.004 + 0.001 

0.39 + 0.02 

4.28 + 0.15 


Table Notes. Col. 1: Normalised fit index; Col. 2 Age; Col. 3: 
Total metal content; Col. 4: Colour excess; Col. 5: Distance from 
the Sun. 


clusters provide some variety of conditions for testing our 
approach. 


4.1 M67 


The high Galactic latitude of M 67 (b = +31.9°) is an advan¬ 
tage for minimising field-star contamination, which makes 
it a convenient cluster for this analysis. Several works have 
come up with a nearly solar meta llicity for M 67 (e.g.|Saraje- 


dini. Potter fc Kirkpatrick|2009 1. The WEBDA[^ database 


provides d© = 0.91 kpc, E{B — V) = 0.06, and the 2.6 Gyrof 
age. More recently, the study of Chen et al.| ( |2014[ ) assumed 
Z = 0.02, (m — M)o = 9.75, and E{B — V) = 0.03, with 
which they derived the significantly older value of 3.5 Gyras 
the the most probable age of M 67. 

Our CMD of M 67 consists of B and V photometry taken 
from the ground-based CCD astrometry catalog of|Yadav et| 


al. 


(20081. We selected stars lying closer than 20' from the 


cluster center, thus resulting in 446 stars brighter than V = 
16. As expected, the evolutionary sequence of a relatively 
old open cluster stands out clearly in the CMD and Hess 
diagram, with a minimum field contamination (Fig. [^. 

Solutions have been searched for isochrones having age 
within 1—13 Gyr and metallicity within 0.001 — 0.03. Our 
approach finds 2 solutions having find > 0.7 for M67. The 
results are given in Table and the corresponding solution 
maps, Hess diagrams and CMDs are shown in Fig. The 
solution map pattern shows a rather tight distribution of 
ages around ~ 3.5 Gyr, but the metallicity distributes over 


® http://www.umvie.ac.at/webda/ 


the relatively wide range 0.06 ^ Z < 0.22. Within the uncer¬ 
tainty, the age we find for M 67 - corresponding to the best fit 
- is nearly the same as that found by Chen et al. (20141, but 
our metallicity is about half theirs. The second-ranked solu¬ 
tion is rather consistent in age, but presents an increase in 
metallicity. Interestingly, the isochrone fit adopted by |Chen] 
et al. polll ) is recovered by our approach as the second- 


ranked solution, with find = 0.82. For comparison purposes. 
Fig. [3] shows the two solutions superimposed to the CMD 
and Hess diagram of M 67. 


4.2 NGC6791 


NGC6791 is closer to the Galactic plane {b = +10.9°) and 
somewhat older and metal richer than M67. For instance, 
WEBDA gives d© = 4.1 kpc, E{B — V) = 0.12, 4.4 Gyr, and 
Z = 0.21 for NGC6791. More recently, [van den Berg et ak] 
(20141 (and references therein) use [Ee/H] ~ 0.3 (thus im¬ 
plying Z ~ 2 Z 0 ) and derive an age of ~ 8 Gyr for NGC 6791. 

The B and V photometry for NGC 6791 was taken from 
the CCD catalogue of [Stetson, Brunt fc Grundahl| ( |2003| ). 
To minimise field-star contamination, we only considered 
stars closer than 3 ' from the center and brighter than V = 
21, resulting in ~ 2100 stars. Even so, the CMD and Hess 
diagram of NGC 6791 clearly contain some degree of field 
contamination (Fig. [^. 

Solutions were searched within the same age and metal¬ 
licity ranges as for M 67. We find 4 solutions having find > 
0.7 for NGC 6791, and they all consistently confirm its older 
age with respect to M67, probably reaching 8 Gyr. Regard¬ 
ing the metallicity, the highest-ranked solutions indicate 
0.019 < Z < 0.028, but with a value as low as Z = 0.012 
with a low probability. Thus, the metallicity appears to be 
higher than solar, probably reaching nearly twice the solar 
value (Table [^. These features (i.e. the relatively narrow 
age range associated with some dispersion in metallicity) are 
present in the solution map. 

The three best-fit PARSEC isochrones produce an ex¬ 
cellent description of the entire evolutionary sequence of 
NGC 6791. However, the difference in find among the first 
three solutions is rather small, which means that they are 
statistically significant as well. The weighted average of the 4 
solutions (using the respective find as weight) yields the age 
7668 ± 193 Myr, Z = 0.021 ± 0.006, E{B - V) = 0.21 ± 0.05, 
and (Iq = 3.90 ± 0.14kpc. An age between 7.5Gyrand 
8 Gyr (in our case derived for most of the main sequence and 
RGB+AGB stars) is consistent with the value obtained by 
Garcfa-Berro et al.|(2010l based on properties of the white 


dwarf cooling sequence of NGC 6791. 


4.3 NGC 2635 


Located in the third Galactic quadrant and at b = +3.96°, 
the field of NGC 2635 is naturally more contaminated than 
both of M67 and NGC 6791. Combined with a large dis¬ 
tance from the Sun and a small population, the available 
CMDs of NGC 2635 lead to somewhat discrepant values of 
its fundamental parameters. For instance, WEBDA provides 
~ 300Myr of age, d© = 5.7kpc, and E{B— V) — 0.35. 
Moitinho et al. (20061 employ CCD broad band photom¬ 
etry and CMDs to find ~ 600 Myr of age d© ~ 4 kpc. 
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E{B — V) = 0.35, and the very-low metallicity Z = 0.004 
for NGC 2635. 

Photometry for NGC 2635 was taken from the CCD 
broad band UBV catalogue of [Moitinho et ^ ( |2006[ ), via 
VizieR. Because of the previous evidence of a younger age 
and lower metallicity (with respect to M 67 and NGC 6791), 
solutions have been searched for isochrones with age within 
lOMyr — 13Gyr and metallicity within 0.001 — 0.02. 

The CMD and Hess diagram of NGC 2635 clearly reflect 
its small stellar content (e.g. the rather discrete distribution 
of stars along the main sequence) and unaccounted-for field- 
star contamination. Nevertheless, our approach indicates ~ 
500 Myr as the most probable age for NGC 2635, together 
with the low metallicity of Z — 0.004. Interestingly, this age 
corresponds to the average between those of WEBDA and 


Moitinho et al. 20061. Metallicity, reddening, and distance 


from the Sun agree with those of Moitinho et al. (20061. The 
best-fit isochrone passes through the main sequence and a 
group of relatively red and bright stars that - in this case 
- would correspond to giants (Fig. [^. It’s obvious that, if 
some field-star contamination is present among the giants, 
the end result of the age may change, but still probably 
within the 100 Myr uncertainty derived by our approach. 
Nevertheless, the find map shows that age and metallicity 
are quite concentrated around ~ 500 Myr and Z = 0.004. 

Finally, we note that the fit index patterns in the so¬ 
lution maps of M67 and NGC 6791 present a high degree 
of similarity. Basically, they are relatively constrained in 
age but somewhat wide in metallicity, which explains the 
^-distribution among the solutions for M67 and NGC 6791 
(Table[^. The apparent degeneracy in metallicity associated 
with isochrones - especially the old ones - is a clear exam¬ 
ple of the difficulties related to deriving precise values of Z 
from photometry. Ideally, isochrones covering a wider range 
of wavelengths - especially towards the blue/violet - should 
mitigate this degeneracy and, as discussed in previous sec¬ 
tions, our approach is designed to work with any isochrone 
set. However, isochrones built with blue filters usually do 
not yet have the same accuracy as those in B and V in 
describing observed evolutionary sequences in CMDs. This 
occurs especially because of the lack of a robust determina¬ 
tion of molecular opacities, incorrect bolometric corrections 
and colour-Te// relations, and incomplete convection de¬ 
scription, mostly related to low-Te//. This issue is explored 
and highlighted, for instance, in Gampos et al. (20131, and 
references therein. 


5 SUMMARY AND CONCLUSIONS 


In this paper we describe a relatively simple and direct ap¬ 
proach designed to extract astrophysical parameters (specifi¬ 
cally the age, metallicity, foreground reddening and distance 
from the Sun) from CMDs of star clusters. Central to the 
approach is the availability of isochrone databases cover¬ 
ing wide ranges - and with a high resolution - in age and 
metallicity. For this reason, we work here with the PARSEC 
( Bressan et al.|[2012 | isochrone database, but any compre¬ 
hensive (in terms of age and metallicity) isochrone set should 
work as well. 

The basic idea is to search for the isochrone that best 
maps the high-density regions of the Hess diagram corre¬ 


sponding to the input CMD. This is accomplished by vary¬ 
ing the apparent distance modulus and reddening for each 
isochrone. At each step, the fit index (find) corresponding 
to the sum of the Hess density spanned by the isochrone 
is computed. The best values of distance modulus and red¬ 
dening are those that lead to the highest possible find for 
each isochrone. Finally, a map of find is produced on the 
ageX metallicity plane, in which the solutions and uncertain¬ 
ties are searched for. 

Tests with simulated CMDs show that the approach is 
efficient in recovering the input values. Obviously, conver¬ 
gence (and uncertainties) depends somewhat on the Hess 
diagram morphology. Star cluster CMDs containing large 
number of member stars tend to produce well-defined evolu¬ 
tionary sequences in Hess diagrams. This, in turn, increases 
convergence of the approach and decreases the uncertainties 
in the retrieved parameters. 

From a theoretical perspective, it’s clear that our ap¬ 
proach depends critically on the accuracy of the isochrones 
in describing the stellar evolutionary sequences at a given 
age and metallicity. Right now this is a minor issue, because 
several groups are actively working on the development of 
isochrones based on state-of-the-art stellar physics. 
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